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ABSTRACT 



Using a sample of 43 bright {V < 16.1, distance < 13 kpc) RR Lyrae stars 

r-| I (RRLS) from the QUEST survey with spectroscopic radial velocities and metal- 

♦p^l licities, we find that several separate halo substructures contribute to the Virgo 

O I overdensity (VOD). While there is little evidence for halo substructure in the 

^ [ spatial distribution of these stars, their distribution in radial velocity reveals two 

moving groups. These results are reinforced when the sample is combined with 

a sample of blue horizontal branch stars that were identified in the SDSS, and 

'i> I the combined sample provides evidence for one additional moving group. These 

^r^ I groups correspond to peaks in the radial velocity distribution of a sample of F 

t^ \ type main-sequence stars that was recently observed in the same directon by 

_ SEGUE, although in one case the RRLS and F star groups may not lie at the 

O ' same distance. One of the new substructures has a very narrow range in metal- 

^-v ' licity, which is more consistent with it being the debris from a destroyed globular 

cluster than from a dwarf galaxy. A small concentration of stars have radial 

rt I velocities that are similar to the Virgo Stellar Stream (VSS) that was identified 

Vh ' previously in a fainter sample of RRLS. Our results suggest that this feature 

extends to distances as short as ~ 12 kpc from its previous detection at ~ 19 

kpc. None of the new groups and only one star in the sample have velocities that 

are consistent with membership in the leading tidal stream from the Sagittarius 

Dwarf Spheroidal Galaxy, which some authors have suggested is the origin of the 

VOD. 
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Subject headings: stars: variables: other 
— Galaxy: structure 



Galaxy: halo — Galaxy: kinematics 



1. INTRODUCTION 

Over the past several decades, a large number of investigations of the stars and the glob- 
ular clusters in the Galactic halo have concluded that the halo formed chaotically and over 
a long period of time through the accretion and disruption of dwarf galaxies (e.g., Searle 
& Zinn 1978, Norris 1986; Carney et al. 1990; Majewski 1992; Mackey & Gilmore 2004; 
see also Freeman & Bland-Hawthorn 2002 for a review). There is evidence, however, that 
the inner and outer parts of the halo have suffici ently different properties that at least two 



mechanisms are requ i red to explain the ha l o fe.g. [Hartwick 



Chiba &: Beers 2000: Kinman et al. 2007: Miceh et al. 



1987 



2007 



Zinn 



199 



Lee et al. 



a 



Maiewski 



2007 



1993 



CaroUo et al. 



20071 ). Accretion appears to have dominated the formation of the outer halo, but the flat- 
tening and prograde rotation of the inner halo and the old ages of the globular clusters in 
the inner halo s uggest that its forma tion involved the dissipative merging of gas clouds at an 
early epoch (see lCaroUo et al.ll2007l . and references therein). While this dual halo picture is 
supported by substantial evidence, it continues to be important to see what fraction of the 
halo can be explained by, for example, an accretion-only scenario. 

The most direct evidence for accretion comes from the detection of halo substructures 
consisting of co-moving streams of stars. The best documented of these are the streams 
originati ng in the Sagittarius (Sgr) dwarf spheroidal (dSph) galaxy, which cross much of 



the sky (lYanny et a. 



Newberg et al. 



2002 



2000 



Ivezic et al. 



2000; Ibata et al. 2001b: Vivas et al. 2001, 2005 



Majewski et al.ll2003l . among others) . Several other substructures have 



been recently detected in the halo of the Milkv Wav (JNewberg et a 



2003 



Belo 



curov et al.ll2006 



2002 



Morrison et al. 



Maiewski et al.l l2003l : lYanny et al.l l2003l : iRocha-Pinto et al. 



among others) and in the halo of M31 



20031 : IZucker et al.ll2004l : iKalirai et al 



200 



i 



2002 



Ibata et al. 



Vivas fc Zinn 



Ibata et al.ll2001a 



2006 



Fardal et al. 



2006 



Ferguson et al. 



2003). These 



streams provide strong evidence in favor of the hierarchical model of the formation of galax- 
ies, in which large galaxies such as the Milky Way formed through multiple mergers with 
smaller building blocks. Simulations of the hierarchical picture p redict that the Milky Way 
has experienced a large number of minor mergers (see for example iBuUock fc Johnstorul2005l ) 
and show that the ones that occurred during the past few Gyrs. may be recognizable to- 
day as stellar streams. While t hese simulations appear to be consistent with the current 
observations of the outer halo (JBell et al.l 120071 ). the search for more streams and better 
documentation of the known ones remains important. 
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We report here an investigation of the halo substructure in the direction of Virgo. 
The first sign of an over de nsity of halo star s in this direction came from the QUEST RR 
Lyrae star (RRLS) survey (IVivas et al.ll200ll ). which revealed a small overdensity near 13 hr 
right ascension (a), which was the west boundary of the survey as it existed then. As the 
survey was extende d, it became clear that this featu re is considerably larger and is centered 
near a ~ 12.4 hr (jVivasI l2002l : I Vivas fc Zinru 120031 ) . In the jargon of the QUEST survey 
it became known as the "12.4 hr dump" (declination (6) was not specified because of the 
narrow range of the survey around 6 = —1°). The RRLS in the clump co yer a range of 
distances from the Sun {D) with the greatest concentration at Z^ ~ 19 kpc (IVivas &: Zinn 
20061 ). Using photometry from the SDSS, Newberg et al. (2002) independently discovered 
an overdensity in F type main-sequence stars in approximately the same direction and at 
a similar distance, which they called S297+63-20. 5. They suggested that it may be related 
to the small group of RRLS identified earlier by I Vivas et al.l (120011 ). and several authors 
have subsequently suggested that the 12.4 hr clump of RRLS and S297+63-20.5 are the 
same substructure or at least related. Also using the SDSS photometric catalogue, Juric 
et al (2005, 2008) later identified a very large overdensity of main sequence stars toward 
Virgo (the Virgo Overdensity (VOD) ). On th e sky, the VOD en compasses both the 12.4 hr 
clump and S297+63-20.5. But, according to IJuric et al.l ((20081), the VOD is spread along 
the line of sight from ~ 6 to > 20 kpc, which is a larger range than eit her the 12.4h clump 
or S297+63-20.5. The spectroscopy of RRLS in the 12.4 hr clump by buffau et all toOdJ ) 
provided the first evidence of a moving group in the Virgo direction, which they called 
the Virgo Stellar Stream (VSS). The r ecent spectroscopy of F type main-sequence stars in 
S297+63-20.5 by iNewberg et al.l (120071 ) revealed a rich structure in velocity space, but with 
a major peak that nearly coincides with the velocity of the VSS. They tentatively conclude 
that S297+63-20.5, the VS S, and the VOD are l i kely associated, if not the same. O n the 
other hand, several authors (JMajewski et al.ll2003l : iMartinez-Delgado et al.ll2C)04l. l2007l) have 



argued that the Virgo substructures are parts of Sgr streams, while others (INewberg et al. 



2OO2I . I2OO7I : IJuric et al.ll2008l : iDuffau et al.ll2006l ) have rejected this hypothesis and suggested 
instead that they are due to a separate merger event involving a dwarf galaxy. 

We are carrying out a program to observe spectroscopically a large number of QUEST 
RRLS in order to measure their radial velocities and metallicities. RRLS are arguably the 
halo tracer that can be most precisely located in space, and with the addition of radial ve- 
locities of even modest precision (±20 km s~ ), they become powerful probes of phase space. 
This wonderful property of RRLS comes with a price, for they are expensive observationally. 
Complete light curves are necessary to obtain both mean magnitudes (hence, precise dis- 
tances) and ephemerides. The later are necessary to isolate the motion of the star in the sky 
from the pulsational velocity of the star by means of the fitting of a radial velocity curve. In 



-4- 



this paper, we present spectroscopic observations of bright RRLS (Vq < 16.1, equivalent to D 
< 12.5 kpc) in the part of the QUEST survey that overlaps with the VOD. After examining 
these data for halo substructure, we discuss the relationships between our results and the 
previous detections of substructures in Virgo and also consider the possibility that the y are 
rel ated to streams from the Sgr dSph galaxy. The RRLS surveys by lKeller et al.l (120071 ) and 
by IWilhelm et al.l (120071 ) provide important data on parts of the VOD that are not covered 
by the QUEST survey. 



2. THE DATA 



2.1. The Sample of RR Lyrae Stars 



According to lJuric et al.l (120081 ) . the VOD occupies ~ 1000 sq. degrees of the sky, and 
its rough limits are shown in Figure [H Notice t hat the SDSS does not extend south of 6 ~ 58° 
at / ~ 300° (see Figure 6 in IJuric et al.ll2008l ). Thus, it is possible that the VOD extends 
farther in that direction. The region surveyed by QUEST is a 2?3-wide band centered at 
S = —1°. It passes over part of the VOD, as shown in the Figure [H The dot-dashed line in 
Figure m shows the southern limit of a second QUEST catalog (Vivas et al. in preparation) 
which is centered at 5 = —3°. Eight stars from this second region are included in our sample. 

The region in common between the VOD and the QUEST survey is roughly limited 
by galactic longitude in the range 260° < / < 340°. At the declination of the QUEST 
survey, those limits are equivalent to 169?5 < a < 212?3. We h ave therefore selected for 
spectroscopy all RRLS in the QUEST catalog (IVivas et al.ll2004j ) within these limits in a, 
and with magnitude V < 16.1. The constraint in magnitude was imposed by the limitations 
of our major instrument, the SMARTS 1.5m telescope (see § 12.21) . The magnitudes of 
this sample correspond to a range of distance from the Sun of ~ 4 to 12.5 kpc (using 
My = 0.55, following Vivas fc Zinn 20 06), which partly overlaps the range estimated for the 
VOD (-6-20 kpc, Ijuric et al.l[2008[ ). 



Figure [2] shows the extinction corrected magnitude Vq (equivalent to distance) of part 
of the QUEST RRLS in the Northern Galactic Hemisphere, as a function of a. The box 
encloses the sample of 44 stars studied in this work. Within the box, all the stars in the 
first QUEST catalog were observed spectroscopically, but only 8 RRLS (out of 38) from the 
second QUEST catalog were observed. For reference, we ha ve indicated with so lid symbols 
the location of the RRLS that probably belong to the VSS (JDuffau et al.ll2006l ). The stars 
in our sample are located in the same area of the sky as the VSS, but are closer to the Sun. 
Most of them (the ones with D > 6 kpc) are located within the same region of the sky and 
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distance range as the VOD. Notice also in Figure [2] that part of the tidal stream of the Sgr 
dSph galaxy lies in the background, at ~ 50 kpc. 



2.2. Spectroscopy 

We obtained a total of 88 spectra of the 44 RRLS with 3 different telescopes between 
2001 and 2007. The instrumental setups of the different observing runs are summarized 
in Table [H and the individual observations are detailed in Table [2l Most of the spectra 
were obtained with the SMARTS 1.5m telescope at Cerro Tololo Interamerican Observatory, 
Chile (denoted as "SMARTS" in Table[2]). We also used observations obtained with the 1.5m 
telescope in La Silla, Chile (denoted as "ESO") and with the Hydra multifiber spectrograph 
at WIYN in Kitt Peak National Observatory, USA. Most of the spectra were taken at blue 
wavelengths, covering the Balmer lines (beginning in H/5) and the Ca II H and K lines. 
However, we took advantage of available time from other projects at the WIYN telescope 
and obtained spectra using a setup which covered the red part of the spectra, including 
the Ca triplet lines. The WIYN observations are denoted in Table [2] as "WIYN-B" and 
"WIYN-R" for blue and red spectra respectively. 

For each target star, we tried to obtain at least two spectra taken at different phases 
during the pulsation cycle. Thirty five out of the 44 RRLS have two or more spectra available. 
Table [2] contains ID (same QUEST ID as in Vivas et al. 2004), Julian date, exposure time, 
telescope, phase of observation (see § 12. 3p . heliocentric radial velocity with its error (see 
§ 12. 5p . and metallicity (see § 12. 4p . The 8 stars from the unpublished second part of the 
QUEST catalog have ID numbers higher than 700. 

The spectrum of star #254 indicates that is too cool to be a RRc star. This is not 
surprising because a small fraction of the QUEST RRc stars are expected to be misidentifi- 
cations of ot her types of variab le stars, mainly variable blue stragglers and W UMa eclipsing 



binaries (see I Vivas et al.ll2004l ). Star # 254 was excluded from any further analysis. 



To avoid excessive broadening in the spectral features due to the changing velocity of the 
star during the pulsational cycle, exposure times were kept short, usually under 30 minutes. 
A few stars have exposure times as long as 60 minutes, which is equivalent to ~ 10% of the 
pulsat ion cycl e of a typical RRLS. In all observing runs, we observed several stars from the 



list of lLaydenI (119941 ) which are both radial velocity and pseudo-equivalent width standards. 
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2.3. Improving the Ephemerides 



Since good eptienierides are essential for fitting tlie radial velocity curves, it is important 
to improve them as much as possible. The QUEST survey has a large numb er of epochs 



(an av erage 32 for our targets), which should yield precise ephemerides. In IVivas et al. 



( 120041 ). the time of maximum light (HJDq) for the RRc stars was set by the observation 
with the brightest magnitude, which is not the best procedure because it can be skewed by 
photometric error and by poor sampling near maximum light. In addition, for some RRc 
stars, it was necessary to list two possible periods because both produced reasonable light 
curve shapes to the eye. 

To obtain better light curve parameters, we have fitted a light curve template to each 
RRc star in our sample. A sine curve produced a poor fit because a typical RRc stars 
has a phase difference between maximum and minimum light of ~ 0.4. Consequently, we 
constructed a light curve template based on 9 of the best observed QUEST RRc stars, which 
did not suffer from period aliases or from a poorly measured HJDq. This template was 
then fitted using ~x^ minimization, while varying the amplitude, phase at maximum light, 
magnitude of maximum light, and period around the values listed in the QUEST catalogue. 
We also explored the parameters around the 1/P ± 1 alias periods. 

The ephemerides for the RRa6 stars in the QUEST catalogue were more secure because 
they were already determined from light curve fitting and because period aliases were less 
of a problem for them. However, the period was not allowed to vary in the fitting. This 
time it was allowed to vary in the fitting, but the resulting light curves are not significantly 
different from the old ones. Only in 5 cases were the changes in period equal to or slightly 
larger than 2 x 10^^ days. The ephemerides that yielded the best fits to the light curves for 
the RRc and RRafo stars have been adopted here and are listed in Table [31 



2.4. Measurements of [Fe/H] 
The metallicities of the R RLS have been measured from the spectrograms following the 



methodology and calibration oflLaydenl (119941 ). which is based on Freeman & Rodgers' (1975) 



modification of the iPrestoru (119591 ) AS method. The method, which involves plotting the 
pseudo-equivalent width of the Ca II K line (W(K)), corrected for interstellar absorption, 
against the mean ps eudo-equivalent w idths of the /?, 7, and b Balmer lines (W(H)), is 



described in detail in I Vivas et al.l (120051 ). 



In general, the type c RRLS are less numerous than the type ab in old stellar populations, 
e.g. globular clusters, and because they are harder to detect than type ab, their number in 
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RRLS surveys of the field suffer more from incompleteness. iLaydenl ( 119941 ) did not produce 
a [Fe/H] calibration for type c RRLS, which were ignored in his study of the metallicities 
and kinematics of field RRLS. There is a sufficient number of type c stars in our sample that 
it is important to measure [Fe/H] from our spectrograms. 

The average type c variable is ~ 600K hotter and has a ~ 0.15 dex larger surface 
gravity than the average type ab variable, although the types overlap in both quantities, 
particularly when the type ab variables are hottest during their pulsations. The Layden 
calibration may therefore require no modification when measuring type c RRLS. This is 
supported by Kemper's (1982) demonstration that the AS method, which was dev i sed fo r 
type ab variables, can be apphed to type c, and more directly by iGratton et al.l (120041 ). 
who found that both types of variables in same globular cluster define a tight sequence in a 
plot of K line strength agains t hydrogen line strength with some intermingling of the types 
along the sequence. Howeve r . iGratton et al.l (120041 ) defined their K line and H line strengths 



differently than did iLaydenl (119941 ). 



To see if Layden's calibration can be applied to the type c, we observed both type c 
and type ab variables in the globular cluster M3 with the Hydra multiobject spectrograph 
on the WIYN telescope. While the fiber setup targeted a large number of variables of 
both types in this variable-rich cluster, the combination of partly cloudy weather and vari- 
able fiber throughputs limited the number of usable spectr a to 8 type c stars and 10 type 
ab. Analyses of high di spersion spectrograms of red giants (ISneden et al.ll2004l ) and RRLS 
( ISandstrom et al.l 120011 ) and the tightness of the red giant branch in the color-magnitude 
diagram of M3 indicate that there is very little if any variation in [Fe/H] among the stars 
in this cluster. For our purposes, we assume that there is no real star to star scatter and 
investigate whether Layden's calibration yields the same [Fe/H] when applied separately to 
the two types of RRLS. 

Our measurements for the M3 RRLS in the Layden pseudo-equivalent width system 
are plotted in Figure [3], where different symbols are used to depict the type ab and type 
c stars. The mean [Fe/H] obtained from all 10 of the type ab variables is —1.78 with a 
standard deviation of 0.27. The large scatter is undoubtedly due to the low S/N of some of 
the spectrograms. The 5 type ab RRLS (filled squares in Figure [3]) for which we have the 
highest S/N spectrograms (50 > S/N > 17 at the Call K line) yield ([Fe/H]) = -1.72 with 
a more reasonable std. dev. of 0.14. The 8 type c variables yield ([Fe/H]) = —1.66 (std. dev. 
= 0.17), and the 5 that have spectrograms with 41 > S/N > 17 yield ([Fe/H]) = —1.63 (std. 
dev. = 0.19). The high S/N samples differ in ([Fe/H]) by 0.09 ± 0.11, which suggests that 
the offset, if any, between the type ab and c variables is small. Since these mean values agree 
well with Zinn & West's (1984) result for M3 ([Fe/H] = -1.66 with an internal precision of 



±0.06), Layden's calibration appears to be equally applicable to type c and type ab RRLS. 

The metallicities that we have derived from the spectrograms of our program stars are 
listed in the last column of Tabled The majority of the spectrograms have S/N > 20, and 
from the experiences with the M3 variables, we estimate that the precision obtained from 
one spectrogram is ±0.15 for the type ab variables and ±0.20 for the type c. The lower 
precision for the type c RRLS is consequence of the convergence of the constant [Fe/H] lines 
at large W(H) (see Figure [3]). A few spectrograms have S/N ~ 15, and the [Fe/H] values 
derived from them are marked with colons in Table [2J Layden's calibration breaks down for 
type ab variables that are observed on the rising branch, where the strength of the K line is 
reduced by the relatively large effective gravity and shock waves may produce emission in the 
cores of the H lines that reduce their equivalent widths. We have discarded measurements 
[Fe/H] that were made at phases where these effects are likely to be important. The surface 
gravity variations of the type c variables are milder, and the equivalent widths of the H lines 
are not affected by emission. None of their measurements have been discarded. 



2.5. Radial Velocities 

Radial velocities for all blue spectra were obtained by Fourier cross-correlation (IRAFs 
fxcor task) with radial velocity standard stars of A and F spectral type, in the spectral 
range 3800 — 5200 A. The radial velocity standards were separated in two groups according 
to their effective temperature. Each target star was cross-correlated with 6 to 14 different 
standard spectra of the group with more similar effective temperature. For each spectrum, 
we obtained a mean radial velocity, V^, by weighting the result from each standard by the 
cross-correlation error returned by fxcor. This error depends on both the S/N of the spectra 
and the similarity between the target star and the standard. 

In order to estimate the radial velocity errors produced by the different telescope- 
spectrograph combinations, we cross-correlated each spectrum of the radial velocity stan- 
dards with the spectra of the other standards, taking care to avoid cross-correlating two 
spectra of the same star. Reassuringly, the differences between the mean velocities from 
these correlations and the literature values were close to zero for all of the standards with 
each instrumental setup. We adopted the standard deviations of these differences as the min- 
imum error in velocity that we can achieve with a particular setup for our target stars (16, 8, 
and 12 km/s for the SMARTS, WIYN blue, and ESO 1.5m spectra, respectively). Because 
these errors include the uncertainties in the zero-points of the wavelength calibrations, we 
believe that these values better reflect the true errors than the means of the cross-correlation 
errors from fxcor , which are smaller. 



-9- 



Cross-correlation was not used to calculate the radial velocities of the red Hydra spectra 
because the subtraction of the strongest night-sky lines left residuals in some wavelength 
regions. Instead, Gaussian line profiles were fitted to the unblended Paschen lines of hydrogen 
and/or the two strongest lines of Ca II, which however at the resolution of the spectrograms 
are each blended with a Paschen line. The Ca II lines were only measured in spectra that 
were taken at the coolest phases of the type ab variables when the Paschen lines are much 
weaker than the Ca II lines. Before computing the radial velocities from the Ca II lines, 
we made slight adjustments to their wavelengths to take into account the blending. These 
were estimated from high S/N spectrograms of bright stars of known velocity whose spectra 
resemble the type ab variables at minimum light. To test this method, we obtained with the 
SMARTS 1.5m telescope red spectra of 6 bright type ab RRLS during cloudy weather. The 
velocities obtained from the measurements of the Ca II lines did not systematically deviated 
from the ones measured from the Paschen lines alone, and the systemic velocities (see below) 
obtained for these stars agreed with values in the literature. For the red Hydra spectra, the 
standard deviation of the mean of the velocities given by the individual lines ranged from 
4 to 10 km s~^, which is consistent with the results obtained from similar spectra of bright 
stars of known velocity. 

The best technique for obtaining the radial velocity of the center of mass (the systemic 
velocity, V^) of a type ab RRLS is to integrate the radial velocity curve given by weak metal 
lines over the whole pu lsational cycle, and the results obtained for RRLS in globular clusters 
(e.g. IStorm et al.lll992l ) are consistent with the radial velocities of the clusters to within the 
internal velocity dispersions of a few km s~^. In our low resolution spectrograms, the weak 
lines produce poorly determined peaks in the cross-correlations, whereas well-defined ones 
are produced by the Balmer lines and other strong features. Furthermore, we could not 
obtain sufficiently large numbers of spectra to determine th e veloci t y curv e for each star. 
For the RRab , we adopted instead the method described by iLaydenl (119941 ). which involves 
measuring the velocities at a few phases (0) from low resolution spectrograms and t hen fitting 
them with the radial velocity curve of the ab variable X Arietis that lOkd (119661 ) obtained 
from measurements of H'j. Observations made near maximum light (0 < 0.1 and (p > 0.85) 
were not used in these fits beca use the radial velocity curve undergoes a large discontinuity 
in this region (see for example ISmithlll995l ). The velocity curve that is determined from 
the H lines has a larger amplitude and di ffers in other ways from the one given by the 



metals lines (see lOke et al.lll962l : IOkdll966l ). Nonetheless, Layden's (1994) comparisons of 



his measurements of Vy with literature values for the same stars indicates that this method 
produces results that are consistent with the errors of the measurements and the fits of the 
velocity curve. In the case of the template star X Ari, the systemic velocity is equal to the 
H line velocity at = 0.5, but there is some small star to star variation in this value. In our 
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calculation of the errors in V^, we have included terms to account for this variation and for 



the variation in the amphtudes of the radial velocity curve (see lVivas et al.ll2005l ). For RRc 
stars, which have much smaller velocity amplitudes than the RRa6 , there are at most small 
offsets between the velocity curves given by the metal lines and the H lines. To measure Vy, 
we fitted a template constr ucted by combining the radial velocity curves of T Sex and DH 



Peg fsee iDuffau et al.l 120061 ) . 



Figure H] shows the best fit for each star having two or more observations. These plots 
include all of the radial velocity measurements for a star, even if they were not used in the 
fits because they are close to = 0. In each case, these points lie within the expected 
range, which provides a check that our measurements do not suffer from large systematic 
errors. Notice also that six of the stars were observed with two different telescopes, and 
the consistency of the measurements, as seen in Figure HI is also an indication that large 
systematic errors are not present. 

In the fitting of the radial velocity curves, the data were weighted by their errors. The 
mean difference between the data points and the radial velocity curve, a^t, in most cases 
was of the same size, or smaller, than the typical errors in the individual measurements (see 
Table [3]). 

For the RRc , we calculated the error in the systemic velocity by combining the individ- 
ual ar in the usual manner for calculating the a of a weighted mean value. For both types 
of stars, if aat was larger than a^ we adopted the former as our final error in the systemic 
velocity. The mean error in the systemic velocities of all, but one, RRLS range from 7 to 
28 km s~^, with a mean value of 17 km s~^. One star, the type c #174, had an unusually 
bad fit of the radial velocity curve and its final error is 42 km s~^. Finally, we calculated the 
radial velocity measured by an observer at the Sun who is at rest with the galactic center, 
Vgsr by 

Vgsr = V^ + 10.0 cos 6 cos / + 225.2 cos 6 sin / + 7.2 sin b; (1) 

Our results are presented in Table [31 which lists: ID, coordinates a and 6, the extinction 
corrected mean V magnitude, the type of RR Lyrae star, the period and Julian date at 
maximum light (see § 12. 3p . distance from the Sun (D), number of spectra available and 
number of spectra actually used in the fitting of the radial velocity curve, the error in the 
fit of the radial velocity curve (erst), the systemic velocity (V^) and its error (a^), and the 
velocity in the galactic rest frame, (Vgsr)- 
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3. DISTRIBUTION OF RADIAL VELOCITIES: SUBSTRUCTURES IN 

VELOCITY SPACE 

As noted above, the VOD encompasses the entire range in a and perhaps a larger range 
of distance than our sample of RRLS. The spatial distribution of our sample of RRLS (see 
the box in Figure [2]) does not reveal any regions of exceptionally high density that are strong 
candidates for halo substructure. There is a roughly uniform band of stars that crosses the 
region at distances between 7 and 10 kpc and also a region of relatively high density near 
12 kpc within 180° < a < 195°. This second feature is the near side of the 



densi t y in RRLS that we have referred to as the "12.4 hr clu mp" (Vivas 



2004 



2002 



arge over- 



Zinn et al. 



Vivas fc Zinnll2006h . which at 19 kpc contains the VSS JPuffau et allbooeh . Fi 



igure 



12 in IVivas &: ZinnI ( 20061 ) shows the location of the 12.4 hr clump in the first band of the 



QUEST survey, but only two of the stars identified there as probable clump members are in 
our sample. The addition of RRLS from the 2nd QUEST band makes the clump even more 
prominent, and there are 9 stars in our sample that now appear to be outliers of the clump. 
If the band of stars at D ~ 8.5 kpc and the 12.4 hr clump are real halo substructures, many 
of the stars should have similar values of Vgsr- 

Figure [5] is a histogram of V^sr of the whole sample of RRLS. The mean velocity is 
+27 km s~^, and the standard deviation is 135 km s~^. These values disagree with recent 
measu rements of samples of BHB stars in the halo (see below). For a sample of 1170 BHB 
stars, ISirko et al.l (l2004b) found a line-of-si ght velocity dispersion, crios = 101.6 km s^^, 
and (V^sr) = km s~^. iBrown et al.l ( 20071 ) found a similar o"ios in their sample of BHB 



selected from the 2MASS survey, although they noted that the velocity dispersion increased 
to 117 km s~^ if they selected only the most metal poor stars i n their sample to avoi d 
thick disk contamination. The Shapiro- Wilk test for normality ( IShapiro fc Wilkl [1965|)lJ 



indicates that the Vgsr distribution of our sample does not deviate significantly from a normal 
distribution. It seems safe, therefore, to use the student t-test and the F-test, which reject 
with > 90% confidenc e that our sample h as the same mean velocity and velocity dispersion 
as the BHB sample of ISirko et al.l (J2004b). However, the velo city dispersion of our sample is 
not significantly different from the value iBrown et al.l (120071 ) obtained for their most metal- 
poor BHB stars. While these differences may be another sign that R RLS and BHB stars in 
the halo have different kinematics (see below and lKinman et al.ll2007l ). the halo substructure 
in our sample is definitely a major contributor to its relatively large velocity dispersion and 
its positive mean velocity. 

While not statistically different from a normal curve, the histogram in Figure [5] is clearly 



^we used the algorithm by iRovstoij (|l995l ). 
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unexpected for a random selection of halo objects, for it has a very prominent peak at Vgsr ~ 
+230 km s~^. In total, there are eight RRLS, or 19% of our sample that have Vgsr > +180 
km s~^, which is very unlikely to be random occurrence (1 in 10,000, if o"ios = 101.6 km s~^ 
and(V^sr) = 0). This excess of RRLS with large positive Vgsr is most likely due to a halo 
substructure. Other likely substructures, which have less extreme values of Vgsr, are revealed 
when distance and velocity are used in concert. 

In Figure El we present histograms of Vgsr for these 3 distance intervals, where we 
have selected dividing lines that avoid putting stars on the boundaries. The histograms in 
Figure [6] do not resemble the normal curves expected of random selections of halo stars, 
and in each case the Shapiro- Wilk test rejects normality with greater than 98% confidence. 
When considering the D < 6.5 kpc interval, it is important to recall that QUEST survey 
is probably incomplete at the bright magnitudes of these relatively near by stars because o f 



saturation of one or more detectors in the QUEST mosaic camera (see IVivas et al.ll2004l ). 
Detector saturation was not a problem for the stars in the more distant intervals, and the 
QUEST survey is expected to be ~ 95% and ~ 70% complete for the type ab and type c 
variables, respectively. For these intervals at least, the departures of the Vgrs distributions 
from the expected halo curves are therefore probably due the presence of halo substructure. 



3.1. Identification of Groups 

To identify possible gr oups of stars, we impleme nted with only a few modifications 



the algorithm described by IClewley fc KinmanI (120061 ) (the Stellar Pair Search). In this 



technique, the distances of each star from all others in the sample are calculated from the 
cartesian Galactic coordinates x, y and z. Pairs of stars are ones that are separated by less 
than a critical distance, dcrtt and have a difference in radial velociticy less than a critical 
velocity, Vcru- Pairs having one star in common form a group. A group may be formed by 
many connected pairs, and the minimum number of stars in a group, A^, is 3 (two pairs 
connected by one star). 

To study the statistical significance of the groups we performed Monte Carlo simulations 
of random samples of stars having the radial distrib ution of the RRLS in the halo. Specif- 



ically, we used the density radial profile described in I Vivas fc ZinnI (120061 ) for the direction 
/ = 300°, b = 60°, which has a slope of —3.1 and takes into account the flattening of the 
inner halo. We assigned to each simulated star a random a and 6 within the same limits of 
our sample, and a random V^sr drawn from a Gaussian distribution with crios = 101.6 km s~ 
and (V^sr) = 0. Each simulated sample was introduced in the group detection algorithm de- 
scribed above. The experiment was repeated many times, and we counted how many times 
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the simulated samples produced groups with characteristics similar to the groups in the real 
sample of RRLS. 

The choices for the parameters dcrit and Vcru require some discussion. As we described 
above, the VOD covers the whole range of right ascension of the sample studied here, which 
is ~ 40°. Any substructure in our data may be as large as this. At 7 kpc from the Sun, 
40° is equivalent to 5 kpc, and at 12 kpc fro m the Sun, it is 8 kpc. On the other hand, the 



uncertainty in distance of the RRLS is ~ 7% (IVivas et al.ll2004j ). which suggests a lower limit 



of dcrit ^ 0.8 kpc. The streams from either destroyed globular clusters or dwarf galaxies are 
expected to be cold substructures, with velocity dispersion < 30 km s~^. Given the errors 
in Vy, Vcrit should be larger than 17 km s~^. Since halo substructures are unlikely to have 
uniform sizes or velocity dispersions, we ran the group detection algorithm and the Monte 
Carlo simulations using combinations of dcrit from 1 to 4 kpc in steps of 0.5 kpc and Vcrit from 
20 to 30 km s^^ in steps of 5 km s^^. A group in the RRLS data was only considered real 
when less than 5% of the 5000 Monte Carlo simulations (a 2a detection) produced a group 
with > A^ members within ±20 km s~^ of |V^sr|, the absolute value of the mean velocity of the 
group of real RRLS. The simulations therefore counted random groups near the Vgsr of the 
suspected group in the real data and ones near —Vgsr- If a group was detected with different 
sets of parameters, as happened in most cases, we report here the parameters that gave the 
lowest probability of random occurrence and/or the maximim number of group members. 
We restricted our sample to the 33 stars that have D > 6.5 kpc since the incompl eteness that 



is cau sed by the different saturation levels of the CCDs of the QUEST camera fIVivas et al. 



20041 ) may effect the interpretation of results at closer distances. 



Two significant groups were found among our RRLS (Figure [7]). As we suspected from 
the histograms in Figures [5] and [6l there is a group with large positive velocity in our data. 
The maximum number of members (8 RRLS) was obtained with dcrit = 3.5 kpc and Vcrit = 25 
km s~^. Not one of the 5000 Monte Carlo simulations produced a similar group. The stars 
making up this group are spread in D from ~ 7.5 to ~ 12.5 kpc (see Figure [7]). Four of 
the more distant stars in this group {D > 11 kpc) are part of the near side of the "12.4 hr 
clump" . The unusually high velocities of these stars and that the fact that some of them are 
probably part of the "12.4 hr clump", unquestionably an over-density in RRLS, suggest that 
they constitute a halo substructure that may be part of the VOD (see §[S]). The combination 
of large Vgsr and a spread in D is consistent w ith a stellar stream that obliquely cuts the 



line-of-sight (see Fig. 15 in lHarding et al.ll200ll ). 



The other significant group contains 9 RRLS. While it is detected with several choices for 
dcrit and Vcrit, the values of 1.5 kpc and 25 km s~^, respectively, yield the lowest probability 
of a random occurrence (only 0.7% of the Monte Carlo simulations produced similar groups. 
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a 2.1a detection). The relatively small value of dcrit of this group indicates that it is more 
concentrated in space than the high velocity group (see Figure [7]). It lies ~ 8.5 kpc from the 
Sun and has a mean Vg^r of —56 km s~^. This group may be due to a stellar stream that 
cuts the line-of-sight at nearly a right angle, which would explain the narrow distribution in 
D and the relatively low VgsT of the stars. 



3.2. Comparison with Blue Horizontal Branch Stars 



If any of the features described above are real substructures, then they should also appear 
in the distributions of the other types of stars that make up very old stellar populations. 
While BHB stars and RRLS are, of course, samples of horizontal branch stars that differ 
primarily in temperature, there are good reasons to suspect that they may not always trace 
the same stellar populations. The majority of the globular clusters in the outer halo have 
red HBs, while the metal-poor globular clusters in the inner halo have primarily blue HBs 
(the dependence of the second parameter of HB morphology on galactocentric distance). 
Consequently, there are large differences i n the frequencies of R RLS and BHB stars between 
the clusters in the inner and outer halos (jSuntzeff et al.lll99ll ). The majority of the dwarf 
spheroidal satellite galaxies of the Milky W ay and M31 have re d HBs and high frequ encies 
of RRLS and low frequencies of BHB stars ( Vivas fc Zinnll2006l ). iKinman et al.l (120071 ) have 
suggested that this variation in the RRLS/BHB star ratio and the kinematic disconnect 
between the inner and outer halos (see § [1]) may be responsible for the kinematic differences 
that they found between samples of RRLS and BHB in the same region of the Galactic halo. 
These observations and the trend in RRLS/BHB suggest that the RRLS may be a denser 
tracer of substructure in the outer halo than BHB stars. Since a small fraction of outer halo 
globular clusters and dSph galaxies have blue HBs, each with minority populations of RRLS, 
and since a small population of BHB stars exists in the red HB clusters, some correlation in 
phase space is expected between these types of stars in the outer halo. On the other hand, 
an exact correspondence would be surprising. 

The sa r nple o f BHB stars considered here was drawn primarily from the survey by 
Sirko et al.l (J2004al ) that is based on photometry and spectroscopy from the SDSS. This 
database of BHB stars is less complete than our RRLS sample because the SDSS spectroscopy 
targeted p rimarily quas a rs an d galaxies. The recent, more complete sample of BHB stars 
stu died bv iBrown et al.l (120071 ) unfortunately does not overlap with the Virgo region. From 



the lSirko et al.l (l2004al ) catalog, we selected all BHB stars (35 in total) within the same range 
of a as the RRLS sample. In 5 we expanded the limits to — 3?5 < 5 < +4?0. This is still a 
small range compared with both the width of the region in a, and the size of the VOD. We 
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also constrained the distance range to 3.0 < D < 13.0 kpc. ISirko et al.l estimate that their 
measurements have an average error in velocity of about 26 km s~^. 

It has been estabhshed both from observations of globular clusters and theoretical mod- 
elling that there is a relation between the luminosity of BHB stars and their te mperature in 



the s ense that bluer BHB stars are fainter than red BHB sta rs (see for example iBrown et al. 



20051 ). When estimating the distances to their BHB sample, ISirko et al.l (l2004al ) adopted in- 
stead the same L/Lq for all BHB stars on the basis of s ome observa tional evidence that they 
realized was in conflict with theoretical computations. ISirko et al.l noted that the adoption 
of the variation in L/Lq with effective temperature that is sugg ested by the comput ations 
would lead to smaller distance moduli by an average of 0.18 mag. iKinman et al.l (120071 ) have 
used the color-magnitude diagrams of globular clusters to derive a relationship between the 
B — V colors of BHB stars and their My values, which is anchored on the same My for 
RRLS that we are using. A com parison, after r nakiri g the necessary transformations be- 
tween photometric systems, of the lKinman et al.l ( 20071 ) relationship to the one employed by 
Sirko et al.l (l2004al ) reveals that the difference in distance modulus is near zero at the blue 
edge of the instability strip and again small at the highest temperatures, but is 0.33 mag (a 
factor of 1.16 in distance) at g — r = —0.23 or B — V = 0. Because only a fraction of the 



sample of BHB stars observed by ISirko et al.l (l2004al ) may be effected by distance errors, we 
have not attempted to recalibrate their distance scale. None of the following comparisons 
depend critically on an exact match of the BHB and RRLS distance scales. 

In a few of the fields that were observed with the Hydra multi-object spectrograph on 
the WIYN telescope, we identified potential BHB stars from the QUEST photometry. The 
three stars that listed in Table H] appear to be true BHB stars. Their {u — g)^ and ((? — r)o 
colors, as measured by the SPSS, are co nsistent with the low gravities of BHB stars on the 
basis of Figure 10 in lYanny et al.l (l2000l ). /^From our spectra, we measured the full widths 
of the /?, 7, and 6 Balmer lines at 0.2 below the continuum level, which are well-know n 
diagnostics of surface gravity at the temperatures of BHB stars (eg. IWilhelm et al.lll999l ). 



In addition to the candidate BHB stars, we measured the WIYN spectra of two BHB stars 
in the globular cluster M3 that are similar in [B — V^)o color to the candidates and of the 
bright A and F type stars that we employed as radial velocity standards with the WIYN 
observations. The candidate BHB stars have smaller H-line widths than the main-sequence 
stars and ones that are similar to the BHB stars in M3. To permit direct comparison with the 
Sirko et al.l (l2004al ) sample of BHB stars, we adopted their method to derive the distances 
to these BHB stars (see Table Hj). 



In Figure El we have plotted histograms of Vgsr for the BHB stars from ISirko et al. 



( l2004al ) and Table H] that lie in the same D intervals as the ones in Figure [61 In this diagram. 
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there may be more smearing between the intervals than for the RRLS because the absolute 
magnitudes of the BHB stars are more uncertain (see above). Neither of the two features 
identified in the RRLS sample is recognizable in the BHB star data alone, which could be 
due to the incompleteness of the sample or to low frequencies of BHB stars in these groups. 
However, several BHB stars have Vgsr values that coincide with the peaks in the RRLS 
histogram. Indeed, when we introduced the combined sample of RRLS and BHB stars into 
the group detection algorithm described above, the same two groups identified in the RRLS 
sample alone appeared. The group with negative velocities now contains 14 members, 5 of 
which are BHB stars. The new mean velocity is —49 km s~^, which we adopt as the velocity 
of the group. Only 5% of the Monte Carlo simulations produced similar groups. The group 
with high positive velocities has 10 members (2 BHB stars). Similar groups were obtained 
in only 2 out of 5000 simulations (0.04%). 

The combined sample of RRLS and BHB stars also contains a third, small but nonethe- 
less significant, group having a mean radial velocity of —171 km s^ . The group, which was 
detected with dcrit = 1 kpc and Vcrit = 25 km s~^, contains only 3 members, all of which are 
BHB stars located in the very narrow intervals oi D = 11.1 to 11.3 kpc and a = 205° to 210°. 
The Monte Carlo simulations found simialr groups in ~ 4% of the cases (a 2.5cr detection). 
This group may be related to main sequence stars with similar velocities (see below) and 
with 2 RRLS at essentially the same distance and velocity but farther west {a ~ 187°). 
Since this group lies near the faint limit of our sample, it may be the near side of a larger 
substructure. 



3.3. Comparison with main-sequence stars 



The recent investigation of the VOD region by lNewberg et al.l (120071 ) has provided radial 
velocity measurements for many stars with the colors of F type stars in two SEGUE fields, at 
(/, b) = (288°, 62°) and (300°, 55°), one of which overlaps with the region we are considering 



here (see Figured]). iNewberg et al.l (120071 ) make the very reasonable assumption that the vast 
majority of these stars are probably near the main-sequence turnoff in the color-magnitude 
diagram. Because such stars have a significant spread in absolute magnitude, their distances 
are much less precisely known than either the RRLS or the BHB stars. The big advantage 
of using these stars is that they are much more numerous than either RRLS or BHB stars. 
Below we examine their data for evidence for the major groups described above. 



seen in our data since very few stars with such high velocities are expected in a sample of 



Group at Vgsr = +215 km s . This group is probably the most remarkable feature 
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halo stars (1 or 2 in our sample of RRLS). The 10 RRLS and BHB stars in this group have 
a mean velocity of {Vgsv) = +215 km s~^, and a standard deviation of 25 km s~^, which is 
not much larger than the typical error in our radial velocities (17 km s~^). As noted above, 
this group appears to be sp read along th e line of sight, and it may be significant that in the 
sample of RRLS studied by lDuffau et al.l (120061 ) . there are 3 stars with Vgsr values consistent 
with membership in this group that have D between 16 and 17 kpc. 



In their discussion of the F-type stars, iNewberg et al.l (120071 ) briefly mention an excess 
of high velocity stars, which w'e suspe ct is related to this group of RRLS and BHB stars. In 
their Figure 9, iNewberg et al.l (120071 ) present for each of their fields histograms of Vgsr for 
the F stars in two ranges in qn ma gnitudes. In th e histo gram for the brighter group, which 
according to the distance scalcl of INewberg et al.l (120071 ) corresponds to 11 to 14 kpc, there 
are 8 F type stars with Vgsr consistent with assignment to our high velocity group. It is 
probably significant that all 8 of these stars are found in the more northern of the two fields 
(see Figure [1]), the one that overlaps with the region studied here (2 stars of even higher 
Vgsr are present in this field as well). Since neither one of the two F-star fields has any stars 
with Vgsr < —200 km s~ , there is an asymmetry toward large positive velocities in the same 
region of the sky where we find one in the RRLS sample. Among the fainter F type stars, 
which lie at distances between 14 and 18 kpc, there are only 3 stars in the 2 fields combined 
with Vgsr within the range of this group. 

We have been able to determine metallicities for 6 of the 8 RRLS in our sample that 
appear to be members of this group. They have a mean ([Fe/H]) = —1.55, with a standard 
deviation of 0.15 dex, which is of the same size as the uncertainty in our [Fe/H] measurements. 
This narrow distribution is much more consistent with the debris from a destr oyed globular 



cluster than from a dwarf galaxy. Two of the 3 high- velocity RRLS measured by lDuffau et al. 



(120061 ) have metallicities that are similar to the metallicities of the brighter RRLS. The other 
one has a much lower ([Fe/H]< —1.95), which may be a sign that it is unrelated. We are 
observing more RRLS in the range 12 < Z^ < 16 kpc, with the expectation that more velocity 
and metallicity measurements will reveal the origin of this substructure. 



^these distance limits and all ones mentio ned in connection to the F-type main-sequence stars have 
substantial unce rtainties. The distan ce scale of INewberg et al.l (|2007f ) is based on the luminosity function of 



Sgr tidal debris (|Newberg et al.ll2002n . which suggests that the assumption of one absolute magnitude for the 
F-type stars yields a distance uncertainty of 16%. If however F-type st ars come from a m ore heterogeneous 
population, resembling for example the globular clusters considered bv lBell et al.l (|2007l ) that span a range 
in [Fe/H], then the distance uncertainty may be as large as 40% 
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Group at V. 



gsr 



—49 km s" 



This is the most populated group that we find in our 



data, containing 9 RRLS and 5 BHB stars with velocities in the range —80 < Vgsr < —10 
km s~^. The mean velocity of this group is (Vgsr) = ~49 km s~^, with a standard deviation 
of 22 km s~^. The contamination from unrelated halo stars is expected to be relatively 
high at the low velocity of this group, and assuming the normal distribution in Figure [5] 
for such field stars, we estimate that only about 50% of the stars in this velocity interval 
are probably members of a stream. Thus, our estimates for both the mean velocity and 
the velocity dispersion are likely to be skewed by this contamination. The 9 RRLS in this 
group have ([Fe/H]) = —1.72, with a standard deviation of 0.28 dex. There is not a strong 
peak in the metallicity distribution that would suggest that the large dispersion is caused 
by contamination by non-members of an otherwise narrow distribution. 



Newberg et al.l (120071 ) found a significant peak at —76 km s in the velocities of the 



F-type stars, which considering the uncertainties is close to the peak that we find. While 
this suggests membership in the same halo substructure, it is not altogether clear that they 
overlap spatiall y. This featur e is rn ost clearly evident in the southern of the two fields 
investigated by iNewberg et al.l (l2007l ). the one that is outside the region investigated here, 
and among the bluer and brighter sample of F-type stars (11 < D < 14 kpc if they have 
Mg = -1-4.2). On the other hand, the RRLS in this group have 7.5 < D < 9.5 kpc, which 
place them in the foreground of the F-type stars, unless the F-type stars are less luminous 
than suspected. There is better spatial coincidence between the F-type stars and a group 



of fainter RRLS that were measured bv lDuffau et al.l (2006), wh ich have similar Vgsr- The 



relationships between the F-type stars, the RRLS in lDuffau et al.l . a s ample of M giants, an d 
the tidal stream from the Sgr dwarf galaxy are discussed at length by lNewberg et al.l (120071 ). 
who suspect that this feature may be a new halo substructure that is unrelated to Sgr. 
The feature that we have identified appears to be unrelated to Sgr as well (se e below). More 



observ ations are required to clarify its relationship to the feature described by lNewberg et al. 



Group at Vr 



gsr 



— 171 km s 



-1 



The 3 BHB stars in the group have basically the same 



values of D, 11.2 kpc, but give n the distance uncertainty a spread of ~ 2 kpc cannot be 



excluded. iNewberg et al.l (120071 ) identify a peak in their F star data at Vgsr = —168 km s 



and in this case, there are good reasons to believe that the features are the same. A peak is 
seen in both of their fields, although it is more prominent in the southern field, the one that 
is outside our region, than in the northern. It is only s een in the brigh t er F stars that have 
0.2 < {g — r)o < 0.4. According to the distance scale of INewberg et al.l (120071 ) these stars lie 
at D between 11 and 14 kpc, which agrees well with the distances of the BHB stars. 
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3.3.1. Does the VSS extends to closer distances? 



The F star data show a very prominent peak at 130 km s"^ in both fields, and lNewberg et al, 
consider this to be the velocity of S297+63-20.5 feature. On the basis of the spatial coinci- 
dence between 8297+63-20.5 and the VSS and the near coincidence in V^..^, they conclude 



that t hese features are probably the same or physically associated. However, iNewberg et al. 



(120071 ) note that the 30 km s" 
explanation. 



offset in Vgsr between S297+63-20.5 and the VSS requires 



The peak at 130 km s~^ i s seen in the faintest a nd the reddest group of F stars, which 
have go between 20 and 20.5. iNewberg et al.l (120071 ) note, however, that the most of these 
stars actually have go < 20.3, which corresponds to D < 16.7 kpc on t heir distance scale. 
Since the core of the VSS is located at about 19 kpc (JDuffau et al.ll2006l ). it is possible that 
30 km s~ offset is due to a gradient in Vgsr with D. Some evid ence for this is provided by 
the data for the brighter and bluer F star data (see Figure 10 in INewberg et al.l ) . which has 
a prominent peak at ~ 150 km s~^. The histograms in Figures El and [H] indicate there are 
5 stars (2 RRLS and 3 BHB stars) having a mean velocity of +143 km s~ which, within 
uncertainties, is identical to the peak seen in the bright F star data. All 5 RRLS and BHB 
stars are located near the far distance limit of our sample at D > 11.8 kpc and within 
the near boundary of the "12.4 hr clump" of RRLS. The 5 stars are indeed identified as a 
group by our algorithm (with parameters dcru = 2.5 kpc, Vcru = 20 km s~^), but by itself 
this group has only low significance because ~ 25% of the Monte Carlo simulations produce 
similar groups. While a chance coincidence cannot be ruled out, the similarity in (Vgsr) 
between this group and the peak in the bright F star data suggests that they are related. 
The group of RRLS and BHB stars provides a hint that the VSS may extend to distances as 
short as ~ 12 kpc. A gradient in Vgsr, with Vgsr becoming less positi ve with increasing -D , 
is expected of a stellar stream that obliquely cuts the line of sight (see lHarding et al.ll200ll ). 
It is also possible, however, that observational error and sample selection may also explain 
the Vgsr offset. We are exploring the possible extension of the VSS and the possibility of a 
gradient by measuring more QUEST RRLS at distances between 13 to 18 kpc (Duffau et al. 
2008, in preparation). 



4. THE RELATIONSHIP TO THE SGR STREAMS 



Soon af ter the first detections of over-densities of stars in th e direction of Virgo, several 
authors (eg. iMajewski et al.l l2003l : iMartinez-Delgado et al.l 120041 ) speculated that they may 
be debris from the disruption of the Sgr dSph galaxy. Even then there was some evidence 
to contrary, for INewberg et al.l (120021 ) had noted that the color of the turnoff in S297+63- 
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20.5 was redder than the color of the turnoff in the Sgr streams. iDuffau et aLl (120061 ) later 
argued that the Sgr streams were not responsible for the VSS or the majority of the other 
stars in their sample because models of the streams produced very few stars at the ob- 



served positio ns and the 



realization bv iJuric et al. 



predic ted values of Vgsr did not match the ones observed. The 



(120081 ) that the VOD has huge dimensions on the sky motivated 



Martinez-Delgado et al.l (120071 ) to reexamine the Sgr connection. Through detailed model- 



ing of the streams, they showed that ones that assume an oblate shape for the gravitational 
potential of the Galaxy can explain the VOD as the leading stream from Sgr as it plunges 
toward the Sun and crosses the galactic plane in the solar neighb orhood. This model does 
not, however, provide a good match to our observations, those of JDuffau et al.l (120061 ) . and 
more generally the 12.4 hr clump of RRLS. 

In the left-hand panels of Figure O we have plotted the RRLS in the first band of the 
QUEST survey that lie within the a range of our sample. Closed sy r nbols depict the RRLS 
that have measured velocities from I Vivas et al.l (120051 ). JDuffau et al.l (120061 ). and the present 
study. These three datasets cover different range of distances, > 40, 16 — 20 and < 1 2 .5 kpc 
respectively. The small crosses are the particles in the stream models of iLaw et al.l (l2005l ) 
under the assumptions that the potential is oblate, spherical, or prolate. The panels on the 
right compare the model predictions for Vgsr as a function of D with the observed velocities. 

The QUEST survey detects the Sgr leading arm at D ~ 50 kpc in the range 195° < a < 
235° (see Figure [2] and Vivas et al. 2005) a nd the 12.4 hr clump in the ranges 12 < D < 20 
kpc and 170° < a < 200° (see Fig. 12 in Ivivas fc ZinnI 120061 ). Note that there is not a 
blending together of these features (see Figure [2]), as one would expect if they are parts of 
the same stream, but instead they are separated by a reg ion of low density . This is not a 
selection effect introduced by the observing technique (see lVivas et al.ll2004l ). 



If the leading arm crosses the galactic plane near the Sun, then the stream stars at small 
D will have large, negative values of Vgsr because they have fallen from ~ 50 kpc on orbits 
that are nearly aligned with the line-of-sight from the Sun. This prediction of large negative 
velocity at small D and the continuit y of the stream are two charac teristics of this picture 
that should not be mode l dependen t . iMartmez-Delgado et al.l (120071 ) note that their models 
reproduce the models of iLaw et al.l (120051 ) . which because they are readily available, we will 
use here to illustrate the salient points of this picture. 



Each of the lLaw et al.l models matches reasonably well the positions and the velocities of 
the RRLS at ~ 50 kpc. However, both the oblate and the spherical models predict smoothly 
varying densities of RRLS over wide ranges in D, which are in conflict with the observe d 
region of relatively low density of RRLS between 20 and 40 kpc (see also lVivas et al.ll2005l ). 
These same models predict that Vgsr becomes more negative as D decreases and reaches quite 
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lar ge \VnKr\ at D < 20 kpc. Only one RRLS in our sample and one in the sample observed 
by lDuffau et al.l (120061 ) have values of Vgsr that are consistent with these models. Also, none 



of the BHB stars in the region have highly negative velocities (see Figure |H]). In the prolate 
model, the leading stream crosses the galactic plane far from the Sun, which explains why 
it predicts relatively few stars in the diagrams at D < 20 kpc. We conclude on the basis of 
these comparisons that none of the halo substructures that we have detected are probably 
due to the leading arm. Moreover, based on the current Sgr models, our observations imply 
that the shape of the Galactic potential can be neither oblate nor spherical. 



Martmez-Delgado et al.l (120071 ) note that in their oblate model the trailing stream of Sgr 
crosses through the region of the VSS and coincides with the VSS in Vgsr (~ 90 km/s). This 
is a plausible explanation for the VSS, if it is still viable after the oblate model is reconciled 
with the observations presented here and with other recent observations whi ch indicate that 



the leading arm does not c ross the galactic plane in the solar neighborhood (JNewberg et al. 
200?! : ISeabroke et al.l[2008[ ). 



DISCUSSION 



Kinematical information is crucial to detect substructures in regions where no obvious 
sp atial overde n sities a re present. Taken together, our RRLS sanaple, th e BHB star sample 
of ISirko et al.l (l2004al ). and the F star sample of iNewberg et al.l (120071 ) paint a consistent 
picture that several separate halo substructures exist at D < 1 3 kpc in the part of the VOD 
that is covered by the QUEST survey. The RRLS sample of iDuffau et al.l (120061 ) and the 
same BHB and F star samples show that some of these substructures extend to larger D. 
None of these substructures seem to be related with the Sgr leading tail. 



The analysis of the SDSS photometric data by iNewberg et al.l (120021 . 120071 ) make a 
compelling case that the region of highest stellar density in S297+63-20.5 lies at Z) ~ 18 
kpc. This is identical to within the err ors of the region of highest density in the 12.4 hr 
clump of QUEST RRLS, which is where iDuffau et al.l (120061 ) discovered the VSS in velocity 



space. The ne a r coin cidence in (Vgsr) between the VSS and the largest peak in Vgsr that 



Newberg et al.l (l2007h observed in S297+6 3-20.5 provides strong evidence that they are the 
same. As noted by INewberg et al.l (120071 ) . the small offset (~ 30 km s~^) could be due to 
measurement error and/or by the small sizes of the samples of spectroscop ically observed 
stars. It is also possible that it is due to a gradient in Vgsr with D (see above). INewberg et al. 



(120071 ) found the same velocity peak in their two spectroscopic fields, which indicates that 
the S297+63-20.5/ VSS has considerable an gular extent o n the sky. Other evidence f or thi s 
is provided by the recent RRLS surveys by iKeller et al.l (120071 ) and IWilhelm et al.l (120071 ). 
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Keller et al.l (120071 ) found 2 clumps of RRLS south of the VOD region (see Figure [T]) at 
coordinates (/, 6) = (301°, 56°) and (324°, 49°) and distan ces 16 and 19 kpc re spectively. 
Velocity information is not yet available for these clumps. IWilhelm et al.l (120071 ) identified 
an excess of candidate RRLS having Vgsr ~ 100 km s~^ in the North Galactic Cap {b > 70°). 
The results presented here indicate that the S297+63-20.5/VSS feature may also extend 
along the line of sight to distances as close as ~ 11.8 kpc. 

We conclude that the VOD contains at least one relatively large structure, S297+63-20.5 
/VSS, which is intersperced with smaller halo substructures that are unrelated to S297+63- 
20.5/VSS or to each other. While these smaller structures might be peculiar to the VOD 
region, it is more likely that they are common in the outer halo. They are better documented 
in the VOD direction simply because the presence of S297+63-20.5/ VSS has drawn attention 
to this region. In the case of our sample of RRLS, the spatial distribution of the stars 
provided only very marginal evidence of any substructure (jVivas fc Zinnll2006l ). The picture 
changed dramatically when we used spatial and velocity data together. There is already 
strong evidence from the spatial distribution of stars in the SDSS t hat much, if not all, of 
the outer halo was accreted from the destroyed of satellite galaxies (jBell et al.l 120071 ). Our 
investigation of the VOD region suggests that the addition of velocity data will make the 
case even stronger. 
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Fig. 1. — Map in Galactic coordinates showing the strip (sohd hnes) surveyed in the first 
QUEST catalog in the Northern Galactic Hemisphere. The dot-d ashed line indicate s the 
region of the second (unpublished) QUEST catalog. According to IJuric et al.l (120081 ) . the 
region of the VOD is roughly enclo sed by the polygon. The small ovals indicate the two 
regions where iNewberg et al.l (120071 ) looked for radial velocity substructures among F stars 
( (lb) = (300°, 5 5°) and (288°, 62°) ). The two x symbols indicate the locations where 



Keller et al.l (120071 ) found overdensities of RRLS at 16 and 19 kpc from the Sun. 
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Fig. 2. — Spatial distribution of the QUEST RRLS in the range of right ascension from 
150° to 240° (from lVivas et al.ll2004j ). The box encloses the bright subsample (Vq < 16.1) 
toward Virgo for which we obtained spectroscopy. The x symbols indicate the 8 RRLS from 
the second QUEST catalog with available spectroscopy that are i ncluded in our anal ysis. 
The large spatial substructures found among the fainter RRLS (see lVivas fc Zinnll2006l ) are 
indicated: the Sgr stream, the globular cluster Pal 5 and its tails, and the "12M clump", 
which contains the VSS. We have indicated with solid circles the RRLS having velocities 
consistent with the VSS (JDuffau et al.ll2006l ). 
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Fig. 3 . — Calibration in the plane W(K) vs W(H) of the metallicity of RRLS, from iLayden 
(Il994l ). Depending on their metallicity, KRab stars lie along one of the dashed lines during 



the pulsation cycle (excluding measurements during the rising branch). Symbols represent 
measurements of RRLS variables in M3, with filled symbols corresponding to high S/N 
observations. Typical error bars for the low and high S/N spectrograms are shown. This 
diagram shows that Layden's calibration is also valid for RRc stars. 
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Fig. 4. — Radial velocity curves for the 33 stars having two or more spectroscopic observa- 
tions. The SMARTS observations are indicated with solid circles, the ESO observation is 
marked with an open triangle, the WIYN-R observations are asterisks, and the WIYN-B are 
open squares. The fitted radial velocity curves are indicated by solid lines. All observations 
are shown including those ones not used in the fitting of the radial velocity curve because 
they were taken near the phase of maximum light (see text). Although they were not used 
in the fitting, they lie within the expected range of velocity. 
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Fig. 5. — Histogram of Vgsr of the 43 RRLS. For reference, the dashed hne indicates the 
expected distribution of velocities in a random halo sample (ahaio = 101.6 km s~^). The size 
of each bin is 30 km s~ . 
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Fig. 6. — Histogram of Vgsr of the RRLS in 3 distance bins. The dashed hne indicates the 
expected distribution of velocities in the halo, normalized to the number of stars in each 
range of distance. 
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Fig. 7. — Distribution of distance from the Sun of the sample of RRLS as a function of right 
ascension. Circles and diamond symbols represent stars with negative and positive radial 
velocities respectively. Colors indicate stars with similar velocities in bins of 50 km s~ . The 
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stars with very high positive velocities and the other one has stars with velocities around 
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Fig. 8. — Histogram of Vgsr of the BHB stars in the same 3 distance bins as Figure [61 The 
dashed hne indicates the expected distribution of velocities in the halo, normalized to the 
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Fig. 9. — Comparison of the locati on and vel o city o f the QUEST RRLS with models of 
the Sgr streams (small crosses) by iLaw et al.l (120051 ) . which assume different shapes for 
the galactic potential (top: oblate; middle: spherical, bottom: prolate). Filled symbols 
correspond to different sets of spectros copic observations of QUEST RR LS: stars with D < 
12.5 kpc (Table 3); 16 < D < 20 kpc jPuffau et al.lboOfih . D > 40 kpc jvivas et al.ll2005l ). 
The open circles in the left panels show the position of QUEST RRLS that have not been 
observed spectroscopically. 



Table 1. Telescopes and Instruments 



Observatory 


Telescope 


Instrument 


Grating 


Spectral Range 


Resolution 


# Spectra/ 


Dates 










(mm^-*^) 


(A) 


(A) 


# Stars 




ESO/La 


Silla 


1.5m 


B&C 


600 


3300- 5500 


3.1 


1/1 


2001 Jun 17 


KPNO 




WIYN 


Hydra 


600 


7100-10000 


2.9 


5/3 


2003 Mar 11 


CTIO 




SMARTS 1.5m 


R-C 


600 


3532-5300 


4.3 


42/28 


2003 May 1-9 


CTIO 




SMARTS 1.5m 


R-C 


600 


3532-5300 


4.3 


24/14 


2004 Feb-Jun^ 


CTIO 




SMARTS 1.5m 


R-C 


600 


3532-5300 


4.3 


5/4 


2006 Jun^ 


CTIO 




SMARTS 1.5m 


R-C 


600 


3532-5300 


4.3 


4/2 


2007 Feb-Apr^ 


KPNO 




WIYN 


Hydra 


400 


3500-6200 


7.1 


9/5 


2007 Feb 20-21 


KPNO 




WIYN 


Hydra 


600 


7100-10000 


2.9 


3/3 


2007 Apr 27, May 2 



^Service mode run 



00 
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Table 2. Individual Observations of the sample of RRLS'' 



ID 


JD 


Texp 


Telescope 


Phase 


Vr 


cr,. 




[Fe/H] 




(d) 


(s) 






(km s'^) 


(km s 


-') 




201 


2452711.7580 


3600 


WIYN-R 


0.76 


58 




5 




201 


2452765.6208 


1200 


SMARTS 


0.02 


44 




16 


-1.95 


203 


2452765.6385 


1200 


SMARTS 


0.02 


13 




16 


-1.40 


212 


2452078.4665 


900 


ESO 


0.56 


88 




11 


-2.11 


212 


2452760.6437 


300 


SMARTS 


0.58 


109 




16 


-1.99 



^The complete version of this table is in the electronic edition of the Journal. 
The printed edition contains only a sample. 

■Uncertain measurements due to low SNR of the spectrogram. This measure- 
ment was not taken into account for calculation of the mean [Fe/H] unless it was 
the only available measurement for the star. 



Table 3. Position, Distance, Velocity and Metallicity of the 43 RRLS*^ 



ID 


a 
(2000.0) 


<5 
(2000.0) 


Vo 


Type 


Period 

(d) 


HJDo 

(d) 


^0 

(kpc) 


N 


Nflt 


(k 


V-, 
ms-i) 


(k, 


o-fit 
ms-i) 


(k, 


tns-1) 


Vgsr 

(km s-1) 


[Fc/H] 


151 


169.625160 


-0.764498 


14.78 


c 


0.29370 


2451611.7002 


7.0 


2 


2 




90 




28 




28 


-35 


-1.32 


157 


171.487530 


-0.161508 


14.39 


ab 


0.70840 


2451571.8154 


5.9 


2 


2 




142 




22 




22 


22 


-1.67 


160 


172.957305 


-2.240574 


15.43 


ab 


0.67619 


2451582.7058 


9.4 


2 


1 




217 








20 


93 


-2.52 


161 


173.844930 


-0.895030 


15.40 


ab 


0.60699 


2451255.6642 


9.3 


2 


2 




18 




7 




15 


-100 


-1.71 



^The complete version of this tabic is in the electronic edition of the Journal. The printed edition contains only a sample. 
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Table 4. Blue Horizontal Branch from WIYN Observations 

ID a S D Vgsr (yvgsr 

(deg) (deg) (kpc) (km s""*^) (km s^^^ 



Q245648 190.865875 -0.926903 9.8 188 8 

Q211556 180.395035 -1.820980 9.3 -6 13 

Q295219 188.952927 -2.765785 12.0 76 8 



